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O . Abstract 

Q^ ' The influence of long-range Coulomb forces on the persistent current of a 
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how to take into account Coulomb effects for the macroscopic quantum tun- 
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q , neling of the Wigner crystal. It is shown that unscreened Coulomb inter- 
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actions increase the stiffness of the Wigner lattice. This leads to additional 

rN I suppression of the amplitude of Aharonov-Bohm oscillations, but makes the 
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predicted anomalous temperature dependence of the persistent current more 
pronounced. The impact of dissipation on the persistent current of a Wigner 
crystal-ring is also studied. 



In condensed matter the Aharonov-Bohm effect manifests itself through magnetic oscilla- 
tions of thermodynamic and kinetic properties of nonsimply connected mesoscopic systems. 
It is a well-known fact that the period of mesoscopic Aharonov-Bohm oscillations is equal 
to the flux quantum $o — hc/e (e is the electron charge) and that the oscillation amplitude 
depends on the properties (conductivity) of the system. 

Until recently the main theoretical and experimental efforts were directed towards inves- 
tigating conductance oscillations in small-sized metallic wires and rings (see e.g. the review 
[1]). Succesful measurements of the flux-induced magnetic moment (persistent current) of 
isolated single metallic- [2] and semiconductor [3] rings have made the study of thermo- 
dynamic properties of mesoscopic quantum rings a central problem of considerable current 
interest. Whereas for metallic rings the problem at hand amounts to properly taking into 
account the effect of impurities on the persistent current of non- (or weakly) interacting elec- 
trons, the main goal in the case of semiconductor rings — formed in a laterally confined high 
mobility two-dimensional electron gas (2DEG) — is to elucidate the role of electron-electron 
correlations in the Aharonov-Bohm oscillations. 

We recall [4] that in a perfect (impurity free) ring even strong Coulomb correlations do 
not influence the persistent current, which is the same as if calculated in the naive model of 
free electrons [5, 6]. This result can be proved as a theorem for Gallilean invariant systems 
(at zero temperature) [7]. In lattice models a drastic change in the magnitude of the persis- 
tent current comes about only in the vicinity of a metal-insulator phase transition induced 
by electron-electron correlations. In the latter case the oscillation amplitude becomes ex- 
ponentially small [8,9,10], which is a characteristic feature of the Aharonov-Bohm effect in 
dielectrics [11, 12]. 

It is reasonable to study the quantum dynamics of electrons in one-dimensional (ID) 
channels -- formed in a 2DEG by a smooth electrostatic potential --in the framework of 
continuum models. In this case Coulomb correlations favor the formation of an ordered 
structure at low density and temperatures. Such a Wigner crystal can exist in rings of large 
enough (mesoscopic) sizes [13]. A simple model of a ID Wigner crystal has been suggested in 



Ref. [14] and has proven to be suitable for studying charge transport in strongly correlated 
electron systems. This model was used in Ref. [15] for an evaluation of the persistent 
current in a ring-shaped Wigner crystal. It was demonstrated that in a perfect ring the 
persistent current caused by the sliding of a Wigner crystal is indistinguishable from the 
current carried by free electrons. The situation is drastically changed in the presence of 
even a small potential barrier. The barrier pins the crystal and charge transport is possible 
only due to the processes of thermal or quantum depinning. At low temperatures the latter 
is dominant and the current becomes of the tunneling type. 

For strong pinning both the zero temperature value and the temperature dependence 
of the oscillation amplitude are significantly changed. An increase in temperature favours 
depinning of the crystal; hence at low temperatures T < T s = hs/L the current increases 
with temperature (s is the sound velocity in the Wigner lattice, L is the circumference of the 
ring). However, in the high temperature region (T > T s ) the persistent current is strongly 
(exponentially) suppressed due a temperature enhancement of destructive interference effects 
(loss of quantum coherence). Thus we get an anomalous temperature behaviour of the 
Aharonov-Bohm oscillations in a pinned Wigner crystal-ring that contradicts the prediction 
of Fermi-liquid theory [6] and is specific for strongly correlated electrons. 

The above mentioned results were obtained in Ref. [15] for a simple model of a Wigner 
crystal without long-range Coulomb interactions and without effects of (possible) dissipative 
motion of the Wigner lattice. The purpose of the present paper is to estimate the influence 
of unscreened Coulomb forces and dissipation on the persistent current of a pinned Wigner 
crystal-ring. 

AID Wigner crystal pinned by a potential V p — local on the scale of the system size L 
— is described by the Lagrangian 

c = S K ~ sV)2 i " Vp5{x){l ~ COSV9) + Z J^ (1) 

Here a is the Wigner lattice spacing (or 1 = n, n being the electron density), m is the 
electron mass, s is the sound velocity (the velocity of charged excitations - plasmons in the 



Wigner crystal), tp = 2nu(x,t)/a stands for the dimensionless displacement field. The last 
term in Eq.(l) describes the Aharonov-Bohm interaction of the Wigner crystal with the 
magnetic flux $. It is worth noticing that the Lagrangian (1) is identical with the bosonized 
form of the Luttinger model (for spinless electrons) if one identifies the parameter r) p , which 
specifies the exponents of different correlation functions in the Luttinger model (see e.g. 
Ref. [16]), with the quantity a = itTi/msa of our model. The parameter a determines the 
"strength" of quantum fluctuations in the system and in what follows it is assumed to be 
small a C 1 (corresponding to a stiff Wigner crystal). 

The long-range part of the Coulomb interaction is decribed by a standard expression, 
which in terms of the displacement field tp looks as follows [17] (see also [18]) 

Ec = £[ Uxdx' d f {x)dMx,) . (2) 



8tt 2 J J ' " yj( x - x'f + d? 

Here d is the length associated with short wavelength "screening" (in fact the width of the 
one-mode channel). We will take it to be of the order of the Wigner crystal spacing d % a. 
In this case plasmons in the Wigner crystal (ka <C 1) experience the influence of unscreened 
Coulomb forces and their dynamics are significantly changed. 

In a constant-potential region the impact of Coulomb interactions, Eq.(2), on the dy- 
namics of plasmons can be accounted for exactly (see e.g. Ref. [17]) 

1 /9 

u(k) = skll + -^—K (ka) I (k ^ 1} s kJ\n(l/ka), (3) 

[ mas 2 J v 

where K (x) is the MacDonald function, Sq = 2e 2 /ma. 

Our prime interest is in the question of how the Coulomb interaction, Eq.(2), influences 
the tunneling dynamics of a pinned Wigner crystal-ring. We will assume that the circumfer- 
ence of the ring is large enough, L ^> a, so that the Coulomb interaction becomes essential 
on scales for which one can neglect the effects of curvature (the distinction between ring-like 
and linear geometry). 

It is physically evident that the main effect of the long-range Coulomb forces is connected 
with the change in the dispersion relation for plasmons, Eq.(3). How will this modification 



influence the tunneling dynamics of a Wigner crystal-ring? It is obvious that Coulomb 
interactions affect the tunneling of the Wigner crystal through the barrier if it is accompanied 
by a distortion of the crystal so that tp' ^ 0. As was shown in Ref. [19], the tunneling 
process is inhomogeneous in the strong pinning regime (in our notations a$Vp ^> T s , where 



T s = hso/L, a — nh/maso — Jclb/cl, clb = h /e 2 m being the Bohr radius). In this case the 
process of depinning the Wigner crystal can be divided into two different stages. In the first 
stage, the field (p rapidly changes its value, Aip = ±27r, in the region l Q ~ amsl/V p <C L that 
contains the barrier. The second stage describes a slow relaxation of the elastic deformation 
that is a result of the first stage; the relaxation is associated with propagation of plasmons. 
It is this second stage that dominates the action in the strong-pinning limit and in our case 
it is modified by the long-range forces. 

The instanton describing the tunneling relaxation is the solution of a free equation of 
motion (in imaginary time r) of the form 

" L/2 . d y (p(y,r) 
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For a screened Coulomb interaction the desired solution is of the form [19] (p = 
±2arctan(sr/|a;|). It is easy to show that for x ~ st ^> a the last term in Eq.(4) can 
be represented in a local form and the modified Larkin-Lee instanton, <Pl-l, obeys the 
equation 
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¥l-l + s ln I — 1 (p L _ L = 0, |x| > a. (5) 

(Recall that we placed the barrier at the point xq = 0). Thus our problem is reduced 
to the one solved by Larkin and Lee [19], but now with an x-dependent velocity s(x) = 



soJ\n(\x\/a). Since s(x) is a smooth function, \s'(x)/s(x)\ ^C 1, we readily get the desired 
instanton solution 



(Pl-l( x , r ) = ±2 arctan | ^ In — - — - (6) 

For the tunneling action corresponding to the trajectory (6) we find (see also [14]) 
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(7) 



The criterion of strong pinning implies that V p ^> e 2 /L. It is reasonable to assume that the 
renormalized barrier height -- renormalized by short wavelength fluctuations — obeys the 



inequality V p < e 2 /a. Then the persistent current at zero temperature takes the form 



Wr = o)~(-i 
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(8) 



where nonessential prefactors have been omitted and N is the total number of electrons in 



the ring. The appearance of the additional factor yhiiV > 1 in the exponent in Eq.(8), can 
easily be understood from the scaling procedure 



s In 1 / 2 ( -) , a' 1 ~ mas -► c^ 1 ln 1 /^^). 



(9) 



Now we consider the temperature behaviour of the persistent current. For a screened 
Coulomb interaction this problem has been studied in Ref. [15]. The trick described above 



(the replacement s — > s(x) = s Jln(\x\/a) ) permits us to use use the results of Ref. [15] 
and readily get the solution. At finite temperature the extremal trajectory is described by 
the function [15] 

7r|x| 
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where (3 = T x is the inverse temperature. The tunneling action corresponding to the 
periodic instanton (10) is 
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where Tj = Hsq/Iq, T c = hso/a. As a function of temperature, the action (11) attains its 
minimum at 



T M ~ 0.5T s Vh[N, A t (Ti 
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(12) 



which leads to the anomalous temperature dependence of the amplitude of the Aharonov- 
Bohm oscillations (see Fig.l). 
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Hence, the unscreened Coulomb interaction results in additional suppression of the per- 
sistent current of a strongly pinned Wigner crystal-ring. The peak in the current plotted as 
a function of temperature therefore becomes more pronounced. These results are consistent 
with the fact that long-range Coulomb interactions, by suppressing the quantum fluctua- 



tions of a Wigner crystal, transform it into an almost classical system a 1 v / ln N — ► oo as 
N ^ oo. 

Now we proceed to the next question. How does dissipation influence the persistent 
current of a Wigner crystal-ring? The first problem is to introduce dissipation in our model 
Lagrangian, Eq.(l). Here we will use a phenomeno logical approach and consider the simplest 
model of dissipative motion, namely linear friction. In other words we will assume that due 
to the presumably weak coupling of electrons in the ring with the reservoir, the motion of a 
Wigner crystal becomes dissipative and can be described in terms of friction in a real time 
formalism. In such an approach it is reasonable to consider the friction as homogenous and 
local. Then for homogenous translations of the Wigner lattice, the dissipative equation of 
motion in real time looks as follows 

mip + rjip — 0. (13) 

Here m is the electron mass and rj is the friction coefficient. This equation, in essence, 
determines the normalization of the phenomenological parameter 77. 

To include friction in the tunneling dynamics, where imaginary time is more appropriate, 
we will follow the standard approach by Caldeira and Legget [20] . According to them linear 
friction in imaginary time, r, can be described by the nonlocal action 

A c ~l ~ V / / drdr 2 , — . (14) 

A) JO ^ sm 2[|( r _ r /)] 

Notice, however, that in our model the trajectories in imaginary time which differ from each 
other by 1-nn (n is an integer) are physically indistinguishable. Therefore the action (14) 
should be slightly modified in order to account for the 27r-symmetry of the model, Eq.(l). 
We propose a dissipative action of the form 
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(15) 
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which leads to the desired form (13) of the equation of motion in real time. It furthermore 
respects the 27r-symmetry of the quantum dynamics of the system in question. It is worth 
noting that the dissipative action (15) is analogous to the one used repeatedly in supercon- 
ductivity for describing the dissipative dynamics of Josephson junctions (see e.g. the review 

[21]). 

For the nonlinear and nonlocal action , Eqs. (1) and (14), we failed to find an exact 
instanton solution. Therefore, in what follows we will restrict ourselves to the case of weak 
dissipation, where a perturbative analysis can be made. It is obvious that the opposite limit, 
the overdamped tunneling "motion" of the Wigner crystal, is described by the trajectories 
with zero winding number (n = 0) and the Aharonov-Bohm oscillations vanish. 

At first we consider the case of a weakly pinned Wigner crystal, aV p <C T s , when the use 
of perturbation theory can be easily justified. For a weakly pinned Wigner crystal which 
tunnels through the impurity homogeneously, it is convenient to treat it as an effective 
particle of mass M = Nm (N is the total number of electrons in the ring) when calculating 
the tunneling action. The characteristic time of tunneling, r*, through the barrier in this 
particular case is much smaller than the periods of time between successive tunneling events, 



t* = a J M/2V P <C (3 [15] and since the effects of dissipation manifest itself most significantly 
on large time scales (~ /3), one can reliably use perturbation theory. 

Substituting unperturbed tunneling trajectories <f(r) = 2-kt/t*, t e [— t*/2,t*/2] into 
Eqs.(l) and (15) one finds 



A t = ^2vr| + 0.3 ^, (16) 

where the first term represents the tunneling action in the dissipationless regime, Tq = 



Hvf/L (V p ^ T ), and the dissipative correction (7 = i]/2m, /17 < T» ~ JV p Tq) can be 
interpreted as being due to the broadening of the quantized levels of the system. 

In the dilute instanton gas approximation (A t ^> 1) the amplitude of the persistent cur- 
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rent is proportional to exp(— At). Hence dissipation exponentially suppresses the Aharonov- 
Bohm oscillations (as it should, see also [22, 12]). Notice that Eq.(16) is valid up to the 



crossover temperature T* ~ JV P T (temperature corrections are of the order T/T*). There- 
fore the criterion for the perturbation expansion to be valid, h'-f <C T*, does not have to be 
violated in the case of moderately weak pinning (V p }t T ^ T* > T ) even if dissipation 
is strong (in terms of free electrons) h'-f •> T . The dissipation, of course, remains weak for 
plasmons hj <T S . 

In the regime of strong pinning, aV p ^> T s , the tunneling of a Wigner crystal is inho- 
mogenous and the characteristic tunneling time is of the order of L/s. In this case one would 
expect a stronger influence of dissipation on tunneling. However, it is easy to show that for 
a stiff (a<l) crystal, the impact of dissipation on the persistent current is similar to the 
case described above. By substituting Eq.(9) into Eq.(15) and doing the integrals one can 
readily get the simple result 

A*-i£ (17) 

Again the perturbative correction due to dissipation does not depend on temperature or 
barrier height. This expression exactly coincides with the dissipative action for an elementary 
transition (A<p = ±27r,n = ±1) of a perfect Wigner crystal 

7i T 9 l h^ . , . 

A > = 2 ¥/ + 2^H («) 

This result can easily be obtained from Eqs.(l) and (15), by making use of the extremal 
trajectories (p n = 2-nnT j 'j3 for a perfect Wigner crystal-ring. One can see from Eq.(18) that 
dissipative corrections are small at temperatures T ^> h'-f. Therefore in the high temperature 
region the impact of dissipation on the persistent current — even for a perfect Wigner crystal 
- can be accounted for perturbatively. The presence of a barrier results in the appearence 
of a new energy scale and this enables us to use perturbation theory even at T — » 0. 

Thus, in all cases considered, the dissipative corrections to the persistent current factorize 
as an overall factor, exp(— const h'-f/To). They depend neither on temperature nor on 
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the strength of the pinning potential (we assume that V p ^ Tivf/L ). This implies that 
dissipation, even if the ring is strongly connected to reservoirs (fry ~ hvp/L), though it 
suppresses the current, does not change the anomalous temperature behaviour predicted in 
Ref. [15]. 
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FIGURE CAPTION 

FIG.l. Temperature dependence of the normalized persistent current in a strongly pinned 
Wigner crystal for two different stiffnesses, a -1 , without long-range Coulomb forces (solid 
line) and when Coulomb forces are considered (dotted line). The Coulomb forces increase 



the effective stiffness of the Wigner lattice (a -1 — ► (Xq 1 VhiN , N = the number of particles. 
Here N = 100) resulting in additional suppression of the persistent current. The peak in the 
current will therefore become more pronounced and the corresponding temperature will be 



shifted by a factor ~ Vln N. Note that the persistent current in all cases is normalized to 
its zero temperature value, which is smaller in the presence of long-range Coulomb forces 

(see text). 
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